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Abstract 

A  fuel  cell  clearly  is  a  multidisciplinary  system.  The  authors  have  chosen  an  original  energy  approach  that  they  have  applied  via  the  Bond  Graph 
formalism.  A  PEM  fuel  cell  Bond  Graph  model  is  proposed  and  detailed.  Its  potentialities  are  then  illustrated  by  studies  on  the  fuel  cell  dynamic 
behavior  (current  harmonic  effects,  etc.)  and  system  simulations  (electricity  generating  unit  based  on  a  PEM  fuel  cell,  etc.). 
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1.  Introduction 

Energy  approach  is  unusual  to  model  electrochemical  com¬ 
ponents  for  both  electrical  engineering  and  electrochemistry 
communities.  Energy  is  however  a  unifying  concept  for  many 
fields  of  physics.  The  electrochemical  components,  especially 
a  fuel  cell,  are  naturally  multi-disciplinary  components  rather 
well-adapted  to  this  approach:  chemistry,  electrochemistry,  ther¬ 
mal  and  electrical  engineering  are  concerned.  The  proposed 
energy  approach  is  applied  via  the  Bond  Graph  modeling. 

The  developed  model  is  an  intermediary  model  between  local 
models  which  solve  differential  equations,  and  models  “black 
boxes”.  The  aim  of  the  authors  is  to  keep  strong  relations  with 
the  physical  phenomena. 

This  energy  approach  has  been  generalized  to  other  electro¬ 
chemical  components  (Li-ion  accumulators,  lead-acid  accumu¬ 
lators,  etc.)  [7]. 

First,  the  authors  recall  the  basis  of  Bond  Graph  modeling. 
Secondly,  they  describe  the  developed  fuel  cell  model.  They 
afterwards  study  the  dynamic  behavior  of  a  PEM  fuel  cell. 
Finally  they  illustrate  some  other  potential  applications  of  this 
model  too. 

2.  Energy  approach 

A  fuel  cell  [3]  clearly  is  a  multidisciplinary  system.  In  order 
to  model  such  a  system,  the  authors  have  chosen  an  energy 
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approach.  This  approach  consists  in  defining  and  modeling  the 
energy  couplings  and  the  energy  exchanges  within  the  system. 
They  have  applied  it  via  the  Bond  Graph  formalism  which  is  an 
explicit  graphical  tool  for  describing  energy  exchanges  within 
a  system.  Allowing  a  unified  representation  of  the  laws  of  the 
various  fields  of  physics,  Bond  Graphs  can  moreover  facilitate 
multidisciplinary  exchanges. 

2.7.  Bond  Graph  principles 

In  Bond  Graphs  [1],  the  energy  exchanges  within  a  system 
are  described  by  bonds  which  represent  the  power  exchanges. 
Two  variables,  effort  and  flow,  are  associated  with  each  bond 
(Fig.  2).  These  factors  have  different  interpretations  in  the  differ¬ 
ent  fields  of  physics  (Table  1).  The  product  of  these  two  variables 
is  the  transferred  power.  The  bond  is  arbitrarily  oriented  by  a  half 
arrow  which  indicates  the  positive  power  flow  orientation.  Fur¬ 
thermore,  causality  is  a  fundamental  concept  in  Bond  Graphs: 
it  defines  the  cause-effect  relations.  The  causal  bar  indicates  the 
effort  direction  (Fig.  2). 

2.2.  Bond  Graph  standard  elements 

Only  a  limited  number  of  elements  are  necessary  to  describe 
the  majority  of  systems.  Table  2  gives  these  main  elements: 
dissipative  element  R ,  inertial  element  I  and  storage  element 
C.  The  connections  between  these  elements  are  implemented 
through  junctions.  There  are  two  types  of  junctions  (Table  3): 
1 -junctions  and  0-junctions  which  respectively  correspond  to 
series  connections  and  parallel  connections.  They  express  in  fact 
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Table  1 

Examples  of  effort  and  flow 


System 

e:  Effort  (unit) 

/:  Flow  (unit) 

Electrical 

v:  Voltage  (V) 

i:  Current  (A) 

Mechanical 

F:  Force  (N) 

V :  Velocity  (ms-1) 

Chemical 

[i:  Chemical  potential 
(J  mol-1) 

dn/dt:  Molar  flow  (mol  s-1) 

Hydraulic 

P:  Pressure  (N) 

dq/dt:  Volume  flow  (m3  s-1) 

Thermal 

T:  Temperature  (K) 

ds/dt:  Entropy  flow  (JK-1  s-1) 

Table  2 

Bond  Graph  elements 

Element 

Represents 

Equation  without  causality 

R:r 

Resistance 

e~rf  =  0 

I:i 

Inertia 

0 

II 

1 

C:c 

Capacitance 

1 

£|^: 

II 

0 

Se 

Effort  source 

e  =  cst 

MS, 

Modulated  effort  source 

e  =  e(input) 

Sf 

Flow  source 

f=  cst 

MSf 

Modulated  flow  source 

/=/(  input) 

the  generalized  Kirchoff’s  laws.  The  transformers  and  the  gyra- 
tors  are  used  to  go  from  a  field  of  physics  to  another. 

Causal  rules  exist  at  the  junctions.  Only  one  port  can  fix  the 
flow  through  a  1 -junction.  Only  one  port  can  fix  the  effort  at  a  0- 
junction.  For  instance,  an  effort  source  connected  to  a  0-junction 
fixes  the  effort  at  the  junction;  the  rest  of  the  system  fixes  the 
flows  through  this  junction. 

3.  Bond  Graph  model  of  a  PEM  fuel  cell 

A  global  scheme  in  Fig.  1  presents  the  different  fields  which 
have  been  modeled,  and  the  main  variables  of  the  model.  Each 


Table  3 

Bond  Graph  junctions 


Junction 

Represents 

Equation 

1 

Equality  of  flows 

y2e‘=° 

0 

Equality  of  effort 

X>=° 

TF 

Transformer 

i 

ei  =  re2,  h  =  rf\ 

GY 

Gyrator 

e\  =  rf2,  e2  =  rf\ 

enerqv  object  e  J  energy  object 

i  — ;  tI  2 

Fig.  2.  Description  of  the  different  fields  of  physics  in  a  PEM  fuel  cell, 
field  of  physics  has  been  represented  by  a  block. 

•  Pgas  represents  the  gas  pressures  (Pn2,  Po2)- 

•  Z)gas  represents  the  gas  volume  flows  ( Z)h2 >  Ad2)- 

•  Et h  is  the  theoretical  thermodynamic  potential. 

•  I  is  the  load  current. 

•  Vpc  is  the  fuel  cell  voltage. 

•  T  is  the  fuel  cell  temperature. 

3.1.  Main  assumptions 

The  main  assumptions  are  as  follows: 


(i)  one-dimensional  modeling; 

(ii)  the  electrodes  are  separately  modeled; 

(iii)  the  input  gases  are  pure  O2  and  pure  H2; 

(iv)  uniform  gas  concentration  in  the  supply  channels  and  no 
pressure  losses  (entry  and  exit  cell  pressures  are  constant 
and  equal); 

(v)  the  gas  diffusion  is  solved  in  steady  state; 

(vi)  no  parasitic  reaction; 

(vii)  the  following  state  parameters  (P,  hydration)  are  known 
and  considered  as  constant  (the  hydration  level  time- 
constant  is  very  low  compared  to  the  studied  phenomenon). 

3.2.  Hydraulic  field 


In  this  field,  the  efforts  are  the  gas  pressures  Pgas  in  bars,  and 
the  flows  are  the  gas  volume  flows  Dgas  (m3  s-1).  The  gas  pres¬ 
sures  are  entry  parameters  of  the  model.  No  losses  are  considered 
in  the  supply  channels.  The  gas  volume  flows  are  calculated  start¬ 
ing  from  the  gas  molar  flow  /;gas  (mol  s- 1 )  and  the  molar  volume 
(m3  mol-1): 


D 


gas 


Tm  J f'gas 


RT 

10  5pgas 


where  R  is  the  perfect  gas  constant. 


(1) 

(2) 
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^Ganode  ( J . ITlOl  )  ^  Fanode  (V)^ 

HTF 

J\  (mol.s  )  /  (A) 

Fig.  3.  Calculation  of  E±  at  the  anode 

3.3.  Chemical  field 


This  part  determines  the  theoretical  thermodynamic  potential 
Et h,  which  represents  the  conversion  of  chemical  energy  to  elec¬ 
trical  energy.  It  is  equal  to  the  difference  between  the  cathodic 
and  anodic  potentials: 


^th  —  ^cathode  ^anode 


(3) 


The  exit  of  the  chemical  block  is  an  effort  (Fig.  4):  the  thermo¬ 
dynamic  potential,  which  results  from  the  free  Gibbs  energy  AG 
at  each  electrode.  The  electrical  circuit  imposes  the  molar  flow 
to  this  chemical  block.  These  transformations  can  be  calculated 
by  the  relations: 


^anode  — 


A  G  anode 

nF  ’ 


^cathode  — 


A  G  cathode 
nF 


/  =  nFfi 


(4) 

(5) 


where  n  is  the  exchanged  electron  mole  number  and  F  is  the  Fara¬ 
day  constant.  The  sign  difference  in  (4)  depends  on  the  reaction 
nature:  an  oxidation  at  the  anode  and  a  reduction  at  the  cathode. 

The  transition  between  the  chemical  field  and  the  electro¬ 
chemical  field  is  realized  thanks  to  a  transformer,  as  shown  in 
Fig.  3.  The  example  of  the  anode  is  given.  The  transformer  ratio 
is  1/nF  as  the  Eqs.  (4)  and  (5)  demonstrate  it. 


3.3.1.  Free  Gibbs  energy  and  dependence  on  gas  pressures 
The  free  Gibbs  energy  variations  are  calculated  from  the  stan¬ 
dard  free  enthalpy  variation  of  each  gas  (H2  for  the  anode,  O2  for 
the  cathode)  at  a  standard  reference  pressure.  They  are  corrected 
by  a  term  taking  into  account  the  pressure  effect.  It  represents 
the  Nernst  law: 

o  RT 

AG  cathode  =  AG  cath0de  ln(7U2)’ 

AGanocie  —  AGano(je  —  RT  ln(PH2)  (6) 


A  ZjO  _  A  ZJ  O  _  A  TjO  ZO\ 

^  ^reaction  —  ^  7  7  products  ^  7  7  reactive  agents 

Enthalpy  and  entropy  variations  of  a  gas  are  temperature  func¬ 
tions: 

A H°(T)  =  A H°(T0)  +  [  ACp  d 0  (9) 

Jt0 

A S°(T)  =  AS°(T0)  +  £  (FF  <W)  (10) 

where: 

=  a  +  p0  +  yd2  (11) 

where  0  is  the  fuel  cell  temperature  and  a ,  /3,  y  are  coefficients 
which  depend  on  the  gas. 

It  is  important  to  note  that  the  power  is  not  conserved 
between  the  entry  and  the  exit  of  this  chemical  block  (Fig.  4): 
PgasT)gas  7^  A//gas/;gas.  Indeed,  the  chemical  energy  transported 
by  the  gas  is  not  taken  into  account  in  the  hydraulic  power.  That 
is  why  all  these  calculations  are  realized  in  the  block  called  “ther- 
mochemical”  in  Fig.  4.  Today,  this  is  not  satisfactory  because 
that  does  not  correspond  to  the  Bond  Graph  spirit.  That  illustrates 
the  difficulty  for  uncoupling  the  thermochemical  phenomena. 


3.3.3.  Produced  reaction  heat 

The  reaction  heat  is  transferred  to  the  thermal  model  as  a 
thermal  power  (“pseudo”  Bond  Graph),  which  is  expressed  in 
Watt  (Fig.  4).  The  temperature  calculated  in  the  thermal  model  is 
injected  into  the  chemical  block,  in  order  to  calculate  free  Gibbs 
energy  variations. 

Note :  The  “pure”  Bond  Graph  modeling  would  handle  tem¬ 
perature  and  entropy  flow.  The  consequence  would  be  that  the 
used  elements  R  would  not  be  homogeneous  to  the  classical  ther¬ 
mal  resistances.  Because  of  this  difficulty,  the  Bond  Graph  users 
handle  classically  the  “pseudo”  Bond  Graph  where  temperature 
is  the  effort  and  the  thermal  power  the  “pseudo”  flow.  The  prod¬ 
uct  “temperature  by  thermal  power”  is  then  not  homogeneous 
to  a  power. 


3.3.2.  Standard  free  Gibbs  energy  variation  calculations 
The  theoretical  recoverable  energy  is  the  enthalpy  variation 
AH  of  the  chemical  reaction.  Unfortunately,  the  practical  recov¬ 
erable  energy,  which  can  be  transformed  into  electrical  energy,  is 
the  standard  free  enthalpy  variation  (standard  free  Gibbs  energy) 
given  by  the  following  relation: 

AG°  =  AH°-TAS°  (7) 

where  A S°  is  the  entropy  variation  (J  mol-1  K-1)  and  TAS°  rep¬ 
resents  the  lost  thermal  energy  (J  mol-1),  which  is  often  called 
reaction  heat. 

The  enthalpy  and  entropy  variations  represent  differences 
between  a  start  state  and  a  final  state  of  the  chemical  reactions. 
They  are  both  calculated  for  each  electrode  thanks  to  this  kind 
of  relation: 


Fig.  4.  Chemical  field;  calculation  of  AG  at  the  anode  (  Jm2  -Ji ). 


R.  Saisset  et  al.  /  Journal  of  Power  Sources  156  (2006)  100-107 


103 


channel  diffusion  layer  reaction  site 


AGano<fe(J.rnOl  ) 

- TF 

J/'(mol.s  ) 


M  A) 


.  Eanode'Ut 


Eanode’Ua 


Double  layer  capacitor^  \g  ^ 


Hi 


TF  ratio  =  1/2F 


*(V) 


Activation  and  Diffusion 
losses 


-H  RS 


Thermal  flow|  T(K) 
(W) 


“J 


Electrolyte  connection 

Thermal  model 

Fig.  6.  Electrochemical  field  at  the  anode. 


3.4.  Electrochemical  field 

The  theoretical  potential  Et h  is  the  fuel  cell  potential  when 
the  load  current  is  zero.  Unfortunately,  this  potential  strongly 
decreases  when  the  fuel  cell  operates.  Different  losses  appear: 
activation,  ohmic  and  diffusion  phenomena.  These  losses  are 
usually  modeled  by  voltage  drops  (overvoltages). 


the  concentration  Cf  reaches  the  value  0: 
( C\A)nFDf 

him  =  “ 


(15) 


The  gas  concentrations  have  a  great  influence  on  the  exchange 
current  of  Butler-Volmer  equation  (12).  The  influence  of  gas 
diffusion  can  be  included  into  this  equation  for  example  at  the 
cathode  [7]: 


3.4.1.  Activation  losses 

At  each  electrode  appear  losses  generated  by  the  reaction 
kinetics  [5].  These  are  called  activation  losses  and  generate 
an  overvoltage  ^act  which  is  subtracted  to  the  theoretical  elec¬ 
trode  potential.  A  non-linear  relation  between  this  overvolt¬ 
age  act,  a  or  c  and  the  electrode  current  /electrode  is  given  by 
the  Butler-Volmer  equation,  if  the  diffusion  phenomena  are 
neglected.  This  law  is  for  example  at  the  anode: 

/anode  =  /o.a  ea^n*F^/RT  -  I0  a  e-(l-«0,a)*a^act,a/«r  (12) 

where  /o,a  is  the  anode  exchange  current  (A),  c^a  is  the  transfer 
coefficient  of  the  oxidation  reaction,  n%  is  the  exchanged  electron 
mole  number  and  iiact,  a  is  the  anodic  overvoltage  referenced  to 
the  anode  potential  (V). 

This  equation  describes  both  the  oxidation  and  reduction 
reactions  that  take  place  at  one  electrode.  Naturally,  a  reaction  is 
predominant  at  each  electrode:  oxidation  at  the  anode  and  reduc¬ 
tion  at  the  cathode.  The  classical  Tafel  equation  can  be  obtained 
with  these  assumptions. 


3.4.2.  Diffusion  losses 

Before  reaching  the  reaction  sites,  the  gas  goes  through  a 
porous  layer,  called  diffusion  layer  (Fig.  5).  The  diffusion  losses 
will  be  higher  at  the  cathode  because  this  electrode  is  the  place 
of  the  water  production.  This  phenomenon  can  be  characterized 
by  the  Fick’s  law,  which  is  solved  here  in  steady  state: 


Df 


d  2c 
dx 2 


(13) 


where  Df  is  the  gas  diffusion  coefficient  through  water  (m2  s  1 ). 
The  final  concentration  at  the  reaction  site  is: 


Cf  =  Ci  - 


Ji8 

~ADf 


—  C\  — 


I  8 
nFA  Df 


(14) 


where  8  is  the  diffusion  layer  width  (m)  and  A  the  cell  area  (m2). 
A  limit  current  I\im  can  be  defined  from  this  expression  when 


exp(((l  -  a0,c)ncF/RT)r]act,c ) 

_  -exp(-(a0,cftc/7^)*?act,c) 

cath0de_  1//0.C  +  l//fim  Ox,c  exp((((l  -  a^)ncF)/RT)rj^c) 
T(l/ Aim  Red,c)  CXp(  —  (fXQ  cncF/ RT)r)&ct,c ) 

(16) 

where  Aim  Ox, c  and  Aim  Red, c  are  the  limit  currents  defined  when 
the  gas  concentration  at  the  reaction  site  reaches  the  value  0. 

3.4.3.  Modeling  of  the  activation  and  diffusion  losses 

The  activation  and  diffusion  losses  are  modeled  by  a  same 
nonlinear  dissipative  element  R  (Fig.  6).  The  relation  which  links 
the  effort  to  the  flow  in  the  Bond  Graph  is  given,  for  example,  by 
(16)  for  the  cathode.  The  associated  losses  represent  a  thermal 
power  which  is  injected  into  the  thermal  model  (Fig.  6). 

3.4.4.  Double  layer  capacitors 

At  each  interface  electrode/electrolyte  appears  the  double 
layer  phenomenon  [2].  Indeed,  no  charge  carrier  transfer  theo¬ 
retically  exists  at  this  interface.  The  distance  between  the  species 
being  very  short,  a  capacitor  of  a  great  value  (several  hundreds 
of  mF  in  our  case)  is  created.  A  constant  element  C  is  used  to 
model  this  phenomenon  (Fig.  6).  This  double  layer  capacitor 
fixes  the  dynamics  of  the  activation  phenomena. 


Electric 

use 


I  I  U anode 

J.  ,, 


Ufuelcell  .  |  U electrolyte  ,  Rl2 (W) 


/ 


/ 


■7  RQl  v' Ar Thermal  model 
no  T  (K) 


U( 


cathode 


Fig.  7.  Connection  of  the  two  electrodes. 
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3.4.5.  Ohmic  losses 

The  ohmic  losses  include  all  the  losses  in  conductors  (elec¬ 
trolyte  and  electrodes)  and  connections.  The  main  losses  being 
a  priori  in  the  electrolyte  (membrane),  the  others  are  not  taken 
into  account  in  our  model. 

A  linear  dissipative  element  R  is  used  to  model  these  losses 
(Fig.  7),  whose  conductivity  depends  on  the  membrane  material. 

For  a  membrane  in  Nation®  1 17,  the  protonic  conductivity  is 
given  by  an  empirical  formula  [4] : 


a  =  (0.005 139A.m  -  0.00326)  e(1267((1/303)“(1/:r)))  (17) 

where  A,m  is  the  hydration  level. 

The  hydration  level  is  determined  experimentally  from  a  mea¬ 
surement  of  the  membrane  resistance  (impedance  spectroscopy 
in  our  case)  and  geometric  data  (membrane  area  and  thickness). 

The  associated  losses  represent  a  thermal  power  which  must 
be  injected  into  the  thermal  model,  like  the  activation  and  diffu¬ 
sion  losses. 


Anode  overvoltage  and 
reaction  heat 


Electrolyte 

losses 


Cathode  overvoltage  and 
reaction  heat 


T 

A 

C  R 


C  R 


ITT 

A;  A 


C  R 


:T 


T  I  . 

1 1  Me 

K  T  !  V  \  V  \  -L 

A  A\AAA 


C  R  1  C  R'c  R'CR 


Extern  sink  1 
Convection  I  Terminal 
and  radiation  i  plate 


Anode 


Electrolyte 

left 


Electrolyte 

right 


Cathode 


I  Terminal 
I  plate 

I 


Extern  sink : 
■  Convection 
land  radiation 


Fig.  9.  Thermal  Bond  Graph  of  an  elementary  PEM  fuel  cell. 
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losses  \  losses  /  losses 


\ 

i  / 

Extern  sink 
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and  radiation 
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plate 

Bipolar 

plate 

Bipolar 

plate 

Terminal 

plate 

Extern  sink 

Convection 
and  radiation 

Fig.  10.  Thermal  model  principle  for  a  three-cell  association. 


The  two  models  of  the  electrodes  are  linked  on  the  level  of  the 
membrane  resistance  (Fig.  7).  The  global  model  is  symmetric 
(Fig.  8). 

3.4.6.  Thermal  field 

Thermal  phenomena  are  very  important  in  a  fuel  cell.  In  a 
first  approach,  two  thermal  flow  directions  are  considered:  lon¬ 
gitudinal  and  transversal.  The  fuel  stack  is  modeled  as  several 
layers.  Each  layer  is  modeled  by  a  RC  element  as  shown  in  Fig.  9. 
The  electrolyte  has  been  divided  into  two  parts  in  order  to  keep 
symmetry,  and  the  losses  are  injected  in  the  middle.  Ambient 
temperature  is  imposed  by  an  effort  source  Se  on  both  sides  of 
the  fuel  cell.  Activation,  diffusion,  and  all  the  ohmic  losses  are 
the  thermal  flows. 

In  a  second  approach,  only  the  bipolar  and  terminal  plates  are 
considered:  the  electrolyte  thickness  (around  200  p,m)  is  much 
thinner  than  the  bipolar  (or  terminal)  plate  (around  1  cm).  The 
thermal  model  is  then  simplified  (Fig.  10). 

The  external  cooling  (fans  in  our  case)  is  taken  into  account 
by  modifying  the  convection  coefficient  (natural  convection: 
5  Wm-2  K-1;  forced  convection:  50  W  m-2  K-1).  The  forced 
convection  coefficient  is  identified  experimentally. 

The  fuel  cell  temperature  is  determined  during  the  model 
simulation. 

4.  Parameter  identification 

The  authors  are  mainly  interested  in  modeling  PEM  fuel 
cell  stacks.  Within  a  stack,  more  or  less  important  disparities 
exist  between  the  different  cells  [6,7].  It  is  then  assumed  that  a 
ft-cell  stack  is  n  times  the  behavior  of  the  equivalent  mean  cell 
[8]. 

The  proposed  model  is  a  model  with  dissociated  electrodes. 
Practically,  and  particularly  in  the  case  of  a  fuel  cell  stack,  only 
the  combined  contribution  of  the  two  electrodes  is  measurable. 
The  authors  have  developed  a  method  to  separate  the  electrodes 
which  will  be  presented  in  a  future  paper. 

The  membrane  resistance  is  first  obtained  from  an  impedance 
spectroscopy  [2,5].  The  impedance  spectroscopy  gives  the 
value  of  the  double  layer  capacitor  too  [7,8].  Electrochemical 
parameters  (activation  and  diffusion)  are  then  identified  starting 
from  a  curve  V(J)  [7,8].  The  thermal  model  parameters  are 
identified  through  thermal  measurements  and  calculations 
(Fig.  11)  [7]. 

5.  Experimental  and  simulation  work 

The  studied  fuel  cell  is  a  commercially  available  200  W  stack 
composed  of  20  cells  (Electrochem  Inc.).  Each  cell  has  an  active 


Fig.  11.  Thermal  measurements  and  simulations  (1  and  3:  natural  cooling;  2: 
fan-cooling). 

area  of  50  cm2.  The  membrane  material  is  Nation®  1 17.  It  is  fan- 
cooled  [8]. 

The  Electrochem  stack  has  been  strongly  exploited  in  our 
studies  [7-9] ,  particularly  about  the  interactions  between  the  fuel 
cells  and  the  static  converters.  Moreover,  the  design  of  a  200  W 
electricity  generating  unit  has  been  studied  by  associating  this 
fuel  cell  with  a  storage  device  [10].  Some  results  are  illustrated 
here. 

The  software  20-sim  was  used  for  all  simulations.  It  is  dedi¬ 
cated  to  the  Bond  Graph  simulation.  Its  use  is  rather  simple  and 
direct. 

5.1.  Dynamic  generation  of  a  voltage-current  curve 

A  curve  V(7)  is  plotted  via  low  frequency  sinusoidal  current 
sweeps.  This  method  is  in  fact  a  dynamic  measurement.  The 
obtained  curve  is  not  the  “static”  curve.  Dynamic  phenomena  are 
thus  solicited  according  to  the  sweep  frequency.  For  example, 
a  “hysteresis”  phenomenon  was  observed  for  low  currents  with 
a  sweeping  frequency  around  1  Hz  (Fig.  12).  The  interactions 
between  the  activation  losses  and  the  double  layer  capacitors 
are  responsible  for  this  hysteresis.  The  behavior  of  the  proposed 
model  is  good. 


Fig.  12.  Experimental  and  simulation  curves  V(l)  in  the  case  of  a  sinusoidal 
current  sweep  (1  Hz). 
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5.2.  Effects  of  high  current  ripples 

A  sinusoidal  current  around  an  operating  point  is  imposed 
to  the  fuel  cell.  In  these  experiments,  the  operating  point  is 
/DC  =  10  A,  the  current  amplitude  Ai  is  3  A  peak  to  peak.  Three 
experimental  curves  are  presented  for  three  different  frequen¬ 
cies:  1  Hz,  100  Hz  and  10  kHz  (Fig.  13). 

A  hysteresis  appears  around  100  Hz  due  to  the  interactions 
between  the  activation  losses  and  the  double  layer  capacitor 
(Fig.  13b). 

This  effect  does  not  occur  around  1  Hz,  because  the  capaci¬ 
tor  has  time  to  charge  and  to  discharge  fully  (Fig.  13a).  In  the 
previous  part,  a  hysteresis  appeared  at  1  Hz,  but  the  range  of  the 


PEM  variable  and 

fuel  cell  controlled 


ultracapacitor 


Fig.  14.  Example  of  a  studied  electricity  generating  unit. 


Fig.  13.  Model  behavior  with  high  current  ripples. 


swept  currents  was  higher  than  in  this  case;  and  particularly  the 
low  currents  are  not  swept  here. 

No  effect  is  observed  around  10  kHz,  because  the  capacitor 
can  be  charged  or  discharged:  the  voltage  across  the  capacitor 
is  constant  (Fig.  13c).  In  other  words,  all  current  harmonics  go 
through  the  double  layer  capacitor  (very  good  filtering).  The 
voltage  ripple  is  only  due  to  the  membrane  resistance.  In  each 
case,  the  Bond  Graph  model  matches  well  the  experimental 
results. 

5.3.  Simulation  application 

The  first  simulation  application  was  the  study  of  a  low  power 
electricity  generating  unit  delivering  a  standardized  single¬ 
phase  voltage  based  on  a  200  W  PEM  fuel  cell  (Fig.  14)  [10]. 
Hybridization  between  a  fuel  cell  and  an  ultracapacitor  was  cho¬ 
sen.  The  idea  was  to  smooth  the  power  delivered  by  the  fuel  cell 
thanks  to  the  device  storage.  Originality  was  the  control  of  the 
fuel  cell  operating  point:  the  fuel  cell  voltage  was  imposed  by  the 
boost  converter  associated  to  the  ultracapacitor  (voltage  pilot¬ 
ing).  Fig.  15  illustrates  the  behavior  of  the  fuel  cell  Bond  Graph 
model  when  a  load  step  occurs.  The  ultracapacitor  supplies  the 
peak  power.  After  the  load  step,  the  fuel  cell  recharges  the  ultra¬ 
capacitor. 

Moreover,  the  authors  have  imagined  other  electrical  piloting 
types  for  a  fuel  cell  (current  piloting,  power  piloting).  They  have 
validated  them  thanks  to  the  fuel  cell  Bond  Graph  model  [7,10]. 
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Fig.  15.  Simulation  of  a  load  step  from  95  to  235  W  during  50  ms — filtered 
powers  and  energy. 
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Finally,  series/parallel  associations  of  fuel  cells  have  been  stud¬ 
ied  thanks  to  this  model  [6,7]. 

6.  Conclusion 

The  energy  approach  does  not  propose  a  local  description  of 
the  physical  phenomena  but  a  global  description  of  their  impact. 
This  model  has  been  realized  in  order  to  be  included  in  electrical 
systems.  The  authors  have  however  demonstrated  in  the  last  part 
that  the  dynamic  behavior  of  a  fuel  cell  can  be  finely  studied  by 
this  approach. 
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